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A method for computing electromagnetic properties of hadrons in lattice QCD is described. The electromagnetic 
field is introduced dynamically, using a noncompact formulation. Employing enhanced electric charges, the 
dependence of the pseudoscalar meson mass on the (anti)quark charges and masses can be accurately calculated. 
At f3 = 5.7 with Wilson action, the tv + — it splitting is found to be 4.9(3) MeV. Using the measured K° — K + 
splitting, we also find m u /md = .512(6). Systematic errors are discussed. Preliminary results for vector meson 
splittings are also presented. 



> 

in 
o 

On 
O 

NO 
ON 

-t— > 
Jh 

Oh! 



X 



1. Light Quark Masses 

If a fundamental theory of quark masses ever 
emerges, it may be as important to resolve the 
theoretical uncertainty in the light quark masses 
as it is to accurately measure the top quark mass 
(e.g. in deciding whether nature avoids the strong 
CP problem via a massless up quark). The par- 
ticle data tables B give wide ranges for the up 
(2 < m u < 8 MeV) and down (5 < m d < 15 
MeV) quarks, while lowest order chiral perturba- 
tion theory gives m u /m d = 0.57 ± 0.04. Nu- 
merical lattice calculations provide, in principle, 
a very precise way of studying the dependence of 
hadron masses on the lagrangian quark mass pa- 
rameters As the electromagnetic contribution 
to hadronic mass splittings within isomultiplets 
is comparable to the up-down quark mass dif- 
ference, an accurate determination of the light 
quark masses requires inclusion of electromag- 
netic effects dressed by nonperturbative QCD dy- 
namics. Here, we discuss a method for studying 
such electromagnetic effects Q. In addition to the 
SU(3) color gauge field, we introduce a U(l) elec- 
tromagnetic field on the lattice which is treated 
by quenched Monte Carlo methods. The result- 
ing SU(3)xU(l) (Coulomb-gauge) configurations 
are then analyzed by standard hadron propagator 
techniques. 



The small size of electromagnetic mass split- 
tings makes their accurate determination by con- 
ventional lattice techniques difficult if the electro- 
magnetic coupling is taken at its physical value. 
We have found that calculations done at larger 
values (roughly 2 to 6 times physical) of the 
quark electric charges lead to accurately measur- 
able isosplittings in the light pseudoscalar me- 
son spectrum, while still allowing perturbative 
extrapolation to physical values. 

We proceed as follows: quark propagators 
are generated in the presence of background 
SU(3)xU(l) fields where the SU(3) component 
represents the usual gluonic gauge degrees of free- 
dom, while the U(l) component incorporates an 
abelian photon field (with a noncompact gauge 
action) which interacts with quarks of specified 
electric charge. Quark propagators are calculated 
for a variety of electric charges and light quark 
mass values. The gauge configurations were gen- 
erated at /3 = 5.7 on a 12 3 x 24 lattice. 200 con- 
figurations each separated by 1000 Monte Carlo 
sweeps were used. In the results reported here, 
we have used four different values of charge given 
by e q =0, -0.4, +0.8, and -1.2 in units in which 
the electron charge is e = ^An/137 = .3028 .... 
For each quark charge we calculate propagators 
for three light quark mass values in order to al- 
low a chiral extrapolation. From the resulting 
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12 quark propagators, 144 quark-antiquark com- 
binations can be formed, leading 78 independent 
meson propagators and masses. 

2. Expected Chiral Behaviour 

Once the full set of meson masses is computed, 
the analysis proceeds by a combination of chi- 
ral and QED perturbation theory. In pure QCD 
it is known that, in the range of masses consid- 
ered here, the square of the pseudoscalar meson 
mass is accurately fit by a linear function of the 
quark masses [||(|. We have found that this lin- 
earity persists even in the presence of electromag- 
netismjlj]. For each of the charge combinations 
studied, the dependence of the squared meson 
mass on the bare quark mass is well described 
by lowest order chiral perturbation theory. Thus 
we write the pseudoscalar mass squared as 

m 2 P = A(e q ,e q ) + m q B(e q ,e q ) + m q B{e q , e q ) (1) 

where e q , eq are the quark and antiquark charges, 
and m q ,mq are the bare quark masses, defined 
in terms of the Wilson hopping parameter by 
— K~ 1 )/2a. (Here a is the lattice spac- 
ing.) Because of the electromagnetic self-energy 
shift, the value of the critical hopping parameter 
must be determined independently for each quark 
charge. This is done by requiring that the mass of 
the neutral pseudoscalar meson vanish at k = k c , 
as discussed below. 

For the physical values of the quark charges, we 
expect that an expansion of the coefficients A and 
B in (|l|) to first order in e 2 should be quite accu- 
rate. For the larger values of QED coupling that 
we use in our numerical investigation, the accu- 
racy of first order perturbation theory is less clear: 
in fact, a good fit to all our data requires small 
but nonzero terms of order e , corresponding to 
two-photon diagrams. Comparison of the order e 
terms with those of order e 2 provides a quantita- 
tive check on the accuracy of QED perturbation 
theory. Only those e terms which significantly 
reduce the \ 2 P er degree of freedom have been 
kept. 

According to Dashen's theorem, in the chiral 
limit the value of m 2 P is proportional to the square 
of the total charge. Thus, we have also allowed 



the values of the critical hopping parameters for 
each of the quark charges to be fit parameters, re- 
quiring that the mass of the neutral mesons van- 
ish in the chiral limit. Thus A takes the form 
A^'(e q + eq) 2 to order e 2 . (Order e 4 terms here 
were not found necessary to fit the data.) The co- 
efficient B in (Q) which parametrizes the slope of 
m 2 P may also be expanded in perturbation theory. 
Of the five possible e 4 terms in B^ 2 \e q , eq), only 
the e 4 , e q eq and e 2 e| terms were found to improve 
the x 2 . The coefficients in A and B, along with 
the four values of k c for the four quark charges, 
constitute a 12-parameter fit to the meson mass 
values. 

3. Lattice Formulation Including EM 

We have chosen a noncompact abelian gauge 
action S cm to ensure that the theory is free in the 
absence of fcrmions, and is always in the noncon- 
fining, massless phase. (Of course, lattice gauge 
invariance still requires a compact gauge-fermion 
coupling). An important aspect of a noncom- 
pact formalism is the necessity for a gauge choice. 
We use QCD lattice configurations which have all 
been converted to Coulomb gauge for previous 
studies of heavy-light mesons. Coulomb gauge 
turns out to be both practically and conceptually 
convenient in the QED sector as well. 

For the electromagnetic action, we take 

Scm = y ^ VjH-^rw — V 'i/A n ^) 2 (2) 

with e the bare electric coupling, n specifies a 
lattice site, the discrete lattice right-gradient 
in the [i direction and A nfl takes on values be- 
tween — oo and +oo. Electromagnetic configura- 
tions were generated using (|J) as a Boltzmann 
weight, subject to the linear Coulomb constraint 
ViA ni = with V a lattice left-gradient opera- 
tor. The action is Gaussian-distributed so it is 
a trivial matter to generate a completely inde- 
pendent set in momentum space, recovering the 
real space Coulomb-gauge configuration by Fast 
Fourier transform. We fix the global gauge free- 
dom remaining after the Coulomb gauge condi- 
tion is imposed by setting the p = mode equal to 
zero for the transverse modes, and the p — mode 
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Table 1 

Coefficients of fitting function, Eq. ([!]). Terms 
consistent with zero were dropped from this fit. 
Numerical values are in GeV 2 and GeV for A and 
B terms respectively. 
Fit 
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0.0143(10)(e 9 + e g -)^ 



1.594(11) 

B^ 0.205(22)e 2 + 0.071(9)e 9 e 5 + 0.050(7)e? 



0.064(17) e 4 - 0.031(4)e 2 e 2 
+0.033(6)e^ - 0.031(4)e 2 e 2 



to zero for the Coulomb modes on each time-slice. 
(This implies a specific treatment of finite vol- 
ume effects which will be discussed below.) The 
resulting Coulomb gauge field A nfl is then pro- 
moted to a compact link variable = e ±lqAni * 
coupled to the quark field in order to describe a 
quark of electric charge ±ge. Quark propagators 
are then computed in the combined SU(3)xU(l) 
gauge field. 

4. Preliminary Results 

For charge zero quarks, propagators were cal- 
culated at hopping parameter 0.161, 0.165, and 
0.1667, corresponding to bare quark masses of 
175, 83, and 53 MeV respectively. The gauge con- 
figurations are generated at (3 — 5.7, and we have 
taken the lattice spacing to be a -1 = 1.15 GeV. 
After shifting by the improved perturbative val- 
ues, we select the same three hopping parameters 
for the nonzero charge quarks [Q. Because this 
shift turns out to be very close to the observed 
shift of k c , the quark masses for nonzero charge 
are nearly the same as those for zero charge. 
For all charge combinations, meson masses were 
extracted by a two-exponential fit (using both 
smeared and local sources) to the pseudoscalar 
propagator over the time range t = 3 to 11. Er- 
rors on each mass value are obtained by a single- 
elimination jackknife. The resulting data is fitted 
to the chiral/QED perturbative formula (Q) by \ 2 
minimization. The fitted parameters are given in 
Table 1. Errors were obtained by performing the 



fit on each jackknifed subensemble. Aside from 
very small corrections of order (m^ — m u ) 2 , the 
7r + — 7T° mass splitting is of purely electromag- 
netic origin, and thus should be directly calcula- 
ble by our method. Because we have used the 
quenched approximation, uu and dd mesons do 
not mix. The squared neutral pion mass is ob- 
tained by averaging the squared masses of the uu 
and dd states. Thus, to zeroth order in e 2 , the 
terms proportional to quark mass cancel in the 



difference 



This difference is essen- 



tially given by the single term 

.2 ~ A (l) e 2 



777?.+ 



(3) 

Using the coefficients listed in Table 1, and 
the experimental values of the ir°,K°, and K + 
masses, we may directly solve the resulting three 
equations for the up, down, and strange masses. 
The 7r + — 7T° splitting may then be calculated, 
including the very small contributions from the 
order e 2 m q terms. We obtain 



4.9 ± 0.3MeV 



(4) 



compared to the experimental value of 4.6 MeV. 
(The electromagnetic contribution to this split- 
ting is estimated @ to be 4.43 ± 0.03 MeV.) Our 
calculation can be compared to the value 4.4 MeV 
(for A QCD = 0.3 GeV and m s = 120 MeV) ob- 
tained by Bardeen, Bijnens and GerardQ using 
large N methods. The values obtained for the 
bare quark masses are 

m u = 3.86(3), m d = 7.54(5), m s = 147(1) (5) 

The errors quoted are statistical only, and are 
computed by a standard jackknife procedure. 
The small statistical errors reflect the accuracy 
of the pseudoscalar mass determinations, and 
should facilitate the future study of systematic er- 
rors (primarily finite volume, continuum extrap- 
olation and quark loop effects) which are ex- 
pected to be considerably larger. The relation- 
ship between lattice bare quark masses and the 
familiar current quark masses in the MS con- 
tinuum regularization is perturbatively calcula- 
ble j3|. For mass ratios (which are independent of 
renormalization prescription) we obtain 



m d - m u 



.0249(3) , — = .512(6) 

m d 



(6) 
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5. Finite Volume Corrections 

The presence of massless, unconfined degrees 
of freedom implies that finite volume effects are 
potentially much larger than in pure QCD, falling 
as inverse powers of the lattice size, instead of ex- 
ponentially. We have estimated the size of these 
corrections phenomenologically along the lines of 
Bardeen, et.al[||, who model the low-g 2 contribu- 
tion to the 7r + — 7T° splitting in terms of 7r, p, and 
Al intermediate states. This analysis gives 

Sm " = 16^J (qi+m{)(q2+m\) dq2 (?) 

We may use this result to estimate the finite vol- 
ume correction by casting the expression as a 
four-dimensional integral over d^q and then con- 
structing a finite volume version of it by replacing 
the integrals with discrete sums (excluding the 
q = mode). For a 12 3 x 24 box with a -1 = 1.15 
GeV, we find that the infinite volume value of 5.1 
MeV is changed to 8m n = 4.8 MeV, indicating 
that the result we have obtained in our lattice 
calculation should be corrected upward by about 
0.3 MeV, or about 6%. In upcoming studies this 
estimate will be checked directly on larger box 
sizes. 

6. Vector Mesons 

The same techniques can be applied to vector 
mesons and heavy-light mesons. For the light vec- 
tor mesons the expected form of the mass matrix 



Table 2 

Coefficients (in lattice units) of fitting function, 
Eq.(|§|) for vector mesons. 

Fit 

~A 0.567(1) + 0.0068(l)(e g + e q ) 2 
B 0.523(47) + 0.205(22) e 2 + 0.073(9)e| 
+0.138(9) e(( e f - 0.027(l) e 2 e | 
-0.017(1)4 - 0.013(1)4 
-0.009(l)4e ? - 0.013(l)e 9 e| 
C 0.81(29) 
D 0.38(50) 



includes the following terms. 

mv = A(e q , e q ) + m q B(e q , e q ) + m q B(e q , e q ) 

+ (m 2 q + m|)C(e g , e q ) + m q m q D(e q , e q ) (8) 

(Including nonanalytic termsj(| (e.g. 0(m 3 / 2 )) 
does not subtantially alter our final results.) Cal- 
culating all the vector mesons mass combinations 
and fitting to this form, we obtain the coefficients 
shown in Table 2. From these results we obtain 
the mass differences 

m p +-m p o = -0.74 ± 1.14(stat) MeV (9) 
m K ,o-m K ,+ = 3.58 ± 0.77(stat) MeV (10) 

7. Conclusions 

Here we have focused mainly on the pseu- 
doscalar meson masses. This is the most precise 
way of determining the quark masses as well as 
providing an important test of the method in the 
7r + — 7T° splitting Q . Further calculations of elec- 
tromagnetic splittings in the vector mesons and 
the baryons Q , as well as in heavy- light systems, 
are possible using the present method. This will 
provide an extensive opportunity to test the pre- 
cision of the method and gain confidence in the 
results. Eventually reliable lattice calculations for 
all isospin breaking effects should be possible. 

We thank Tao Han, George Hockney, Paul 
Mackenzie and Tetsuya Onogi for contributions 
to our effort. 
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